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The properties of vortices in superconducting thin films are revisited. The interaction between 
two Pearl vortices in an infinite film is approximated at all distances by a simple expression. The 
interaction of a vortex with a regular lattice of real or image vortices is given. The two spring 
constants are calculated that one vortex in the vortex lattice feels when the surrounding vortices 
are rigidly pinned or are free. The modification of these London results by the finite size of real 
films is obtained. In finite films, the interaction force between two vortices is not a central force 
but depends on both vortex positions, not only on their distance. At the film edges the interaction 
energy is zero and the interaction force is peaked. Even far from the edges the vortex interaction 
considerably deviates from the Pearl result and is always smaller than it. 

PACS numbers: 74.25.Qt, 74.78.Db, 74.78. Bz, 74.25. Ha 



; I. INTRODUCTION 

In this note the problem of vortices in thin supercon- 
, ducting films is revisited. Useful formulae for the in- 

■ teraction between vortices in thin superconducting films 
within London Theory are given when the film thick- 

■ ness d is smaller than the London penetration depth 
[ A. In such thin films the vortex interaction is mediated 
. mainly by the magnetic stray field, and the screening 
' of the in-plane supercurrents is governed by the effec- 
, five 2-dimensional (2D) penetration depth A = X^/d > 
' X > d. The current density j{x,y,z) in thin films is 

nearly independent of the z coordinate (perpendicular 
to the film) and the sheet current in the film is thus 

3{x,y) — /^^/2 2/' '^)'^'^ ~ id. The sheet current of a 

I vortex in such a thin film and the interaction force be- 

. tween two vortices was calculated first by Judea Pearl^ 

■ for an infinitely large film. For infinite films of arbi- 
. trary thickness the Pearl-London vortex and the vortex- 
' vortex interaction were calculated analyticallyj^i^!^ and 
, the Ginzburg-Landau (GL) theory for periodic vortex lat- 
' tices in such films with arbitrary GL parameter k in the 
' entire range of inductions B was computed in Ref. [^. 

] For superconductors of finite size, recently there have 

■ been numerous computations based on the static or time- 
dependent GL equations for small (mesoscopic) speci- 
mens containing a few vortices, also giant vortices with 
several flux quanta, and considering the equilibrium state 
or the penetration, exit, nucleation, and annihilation of 
vortices. For examp le, sup erconducting disks are consid- 
ered in RefJ6|j7p^[l3^^ squares and other shapes 
in Ref. [T3lll4lll5lll6l . also squares containing magnetic 
dotsiiii^ or antidots or blind holes^SiiS as pinning centers, 
and disks in inhomogcncous magnetic field. When the 
2D penetration depth A = A^/d is much larger than the 
film size, the energy of the magnetic stray field outside 
the specimen is expected to be negligible and analytical 
approximations can be madci^ Further work computes 
the properties of infinitely long straight vortices in semi- 
infi.nitc^'^ or cylindrical2i samples and argues that these 



results apply also to thin films in a perpendicular field 
in cases when the magnetic stray-field energy outside the 
film may be disregarded. 

The present paper focusses on situations when the 
stray-field energy is not negligible, e.g., since A is smaller 
than the film size, and it is restricted to the London limit 
in which the superconducting coherence length ^ and vor- 
tex core radius Tc « ^ are much smaller than the film size 
and than A. In Sec. II we first consider the limit of ideal 
screening (A = 0) when the sheet current and energy of 
vortices are completely determined by the stray field and 
can be derived easily for an infinite film. Introduction of 
a finite A then leads to Pearl vortices, whose interaction is 
discussed in some detail in Sec. II-V and in Appendix A. 
Finally, in Sec. VI vortices and their interaction in films 
of finite size are presented, computed by a metho d^^i^^ 
that applies to thin flat films of any size and shape and to 
any A. As a main result in Sec. VII it will be shown that 
the sheet current and interaction of vortices in films of 
finite size differ considerably from the ideal Pearl result. 
As expected, this deviation is most pronounced near the 
film edges, where the real interaction potential vanishes 
while the interaction /orce is peaked. But less expectedly, 
this deviation is considerable even when the interacting 
vortex pair is located far from the film edges, e.g., near 
the center of a square film: The real interaction with a 
central vortex is reduced from its Pearl value by a fac- 
tor which decreases approximately linearly from unity to 
zero as one goes from the central vortex towards the edge, 
see Fig. 6 below. A summary is given in Sec. VIII. 



II. VORTICES IN THIN INFINITE FILMS 

The thin-film problem differs from the behavior of cur- 
rents and vortices in bulk superconductors by the dom- 
inating role of the magnetic stray field outside the film. 
The interaction between vortices occurs mainly by this 
stray field and thus has very long range (interacting over 
the entire film width) , while in bulk superconductors the 
vortex currents and the vortex interaction are screened 
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and thus decrease exponentiaUy over the length A. 

The long-range currents and forces can be shown for 
the ideal screening limit of zero A as follows. Consider 
one vortex in the center of a large circular film with ra- 
dius R oo. On length scales larger than A this point 
vortex behaves like a magnetic dipole, composed of two 
magnetic monopoles: one sends magnetic flux $o into the 
space above the film and one receives the same flux from 
the lower half space. Here $o = h/2e = 2.07-10-15 Tm^ 
is the quantum of magnetic flux. The magnetic field lines 
of this point-vortex are straight radial lines, all passing 
through this point. The magnitude of this magnetic stray 
field is $o/(27rr|) above and — $o/(27rr§) below the film, 
since the flux <i>o has to pass through the shell of any half 
sphere of 3D radius r3 — {r^ + z^Y^^ and surface area 
27rr|. Here r — {x^ + y^)^/^ is the 2D radius in the fllm 
plane z — 0. Directly above and below the fllm plane 
the stray fleld has the same magnitude ^o/{2Trr'^) but 
has opposite sign. This jump of the magnetic fleld com- 
ponent parallel to the fllm is caused by a sheet current 
that circulates around the vortex and equals this fleld 
difference in size, J(r) = ^o/diQirr'^). This current thus 
decreases very slowly with radius r. The force F between 
this central vortex and a second vortex at a distance r 
is radial (central force) and repulsive and in size equals 
just this sheet current J{r) times the flux quantum ^q, 



F(r) = $0 J(0 



$2 



(1) 



The interaction potential V{r) between these two point 
vortices (Pearl vortices) follows from its derivative 
F{r) = —V{ry and the deflnition V{oo) = 0, 



V{r) 



(2) 



Note that this result does not depend on A = X^/d. It 
applies also for finite A if the distance is large, r ^ A. 

The general results for arbitrary A can be derived from 
the expression for the interaction potential noting that 
stifl F{r) = -V{ry and J(r) = F{r)/^o- One has 



V{r) = 



d^k 2 cos kr 
47r2 k + 2Afc2 
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'dk 2Ja{kr) 



27r 1 + 2Afc 



(3) 



fc2. The 



flO J 47r^ + Z1\H^ f^O JO 

where Jo (a;) is a Bessel function and fc2 = k'^ 
limits for small and large distances are,— see also Eq. (2) 

ln(2.27A/r 



Vir) 



27rA 



for r < A , (4) 
for r > A . (5) 



The factor 2.27 is a fitted constant that follows from the 
numerical evaluation of the integral (3). An excellent fit 
to the exact result valid for all distances < r < oo and 
all A is within line thickness 
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Figure 1 shows the exact potential from Eq. (3) , its limits 
(4) and (5), and the approximation (6). One can see that 
the expression (6) (dots) practically coincides with the 
exact potential (solid line). 

In the original paper by Pearl^ the force fi2{'r) — 
— V'(r) between two point vortices was expressed in 
terms of Struve and Neumann functions Si (or H„) and 
Ni (or Yn, the Bessel function of the second kind, or 
Weber function))^ 



fi2{r) 



$2 



8A2^o 



«(s) 
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Though both 5*1 and A'^i are oscillating functions this 
force /i2 is monotonic and agrees with —V'{r), Eq. (3). 
The interaction potential V{r) may be also expressed in 
terms of Struve and Weber functions, see appendix in^S,, 
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III. INTERACTION WITH VORTEX LATTICES 

Next I consider the interaction potential Vir) between 
a vortex at position r = (x, y) in the film and a regular 
lattice of vortices sitting at ideal-lattice points R with 
reciprocal lattice vectors K. This situation can be inter- 
preted in two ways: 

1. This lattice can consist of real vortices in the film, 
generated by a constant applied magnetic field Ba per- 
pendicular to the film plane (along z). In this case it will 
be typically a triangular lattice with density n = Ba/^o, 
since the magnetic field has to penetrate the infinitely 
extended film completely, i.e., the average induction B 
equals Ba- The lattice spacing is then a ~ (<i>o/i?)i/^. 

2. In a second application the lattice may be formed 
by the image vortices sitting on a superlattice that is con- 
structed to achieve periodic boundary conditions in the 
numerical simulation of a film with vortex pinning. In 
this case the superlattice may be chosen such that the 
basic cell is a rectangle with side lengths iLx,Ly). The 
regular lattice positions are then R — ifiL^, vLy) where 
(1,1^ = 0, ±1, ±2, . . . are integers, and the reciprocal lat- 
tice vectors are K = i2n ^ / Lx,2'kv / Ly). 

The interaction potential between the probing vortex 
at position r and the periodic lattice of vortices (or im- 
ages) is obtained by linear superposition of potentials of 
the form (3), 



^pcr(r) 



d2fc ^ 2exp[ik(r- R)] 



47r2 ^ 

R 



k + 2Afc2 



(9) 



This potential is periodic, Vpor(r-f R) — Vpcr(r). The 
infinite lattice sum in (8) can be evaluated using the for- 
mula 



E 

R 



exp(ikR) 



47r^n ^ 82 

K 



k-K) 



(10) 
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FIG. 1: The interaction potential V{r), Eq. (3), between two 
Pearl vortices (solid lines), the limits (4) and (5) (dashed 
lines) and the approximation (6) (dots). V is in units 
$Q/(/ioA) and r in units A = X^/d. The upper plot shows 
directly V{r). The lower plot shows VV versus ^/r to de- 
pict larger ranges of V and r and to show the validity of the 
limiting expressions at small and large r more clearly. 



where n is the density of lattice points R; the density of 
reciprocal lattice points is l/(47r^ri). Inserting (10) into 
(9) and performing the k integration over the 2D delta 
functions (52 (k — K) one obtains 



V^pcr(r) = 



$0 n 
Mo 



2 cos(Kr) 



(11) 



This periodic potential may be used for computer simula- 
tion of vortex pinning in a thin film with periodic bound- 
ary conditions. In this case the super- lattice density is 
n = l/{LxLy). 

Next I apply this potential to the following problem. 
Consider a regular (e.g. triangular) vortex lattice in the 
film in which all vortices are pinned to these ideal vortex 
positions and the central vortex at R = is moved a 
bit in X direction. Which interaction potential with all 
other vortices does this central vortex see? To obtain this 
potential I4(r), or its curvature V^'{Qi) = d'^Vcir) / dx^ at 



r — 0, we have to subtract from the complete lattice 
sum (9) the term with R = 0, since this term describes 
the interaction of the shifted vortex at r « with the 
unshifted vortex at R = 0, which would be very large 
and diverging as r — s- R. The resulting potential is 



R#0 



Mo 



2 cos(Kr) 
K + 2AA'2 



2exp[ik(r - R)] 
k + 2Afc2 

[ cPk 2cos(kr) 



47r2n k + 2Afc2 
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with B ~ n^Q. This result is still general and applies 
to any r. To obtain its curvature at r = we apply to 
it the Laplace operator = d"^ jdx^ -I- d"^ jdy^. This 
creates a factor in the sum and a factor — fc^ in 

the integral. Due to symmetry of the triangular (and 
also the square) vortex lattice, the curvature of V^(r) 
at r = is the same along all directions, i.e., one has 
K(r) = K(0) + iT/"(0)r2-HO(a;4,a;V,y') with ^"(0) = 
d'^Vcldx'^ = d'^Vjdy'^ = iV^y, at x = y = 0. The resuh 
for this curvature is thus 



K"(0) = 



£oS 
2/^0 



2fc2 



(fk 



A-K^n k+ 2Ak^ 



K#0 



2K^ 



K+2AK^ 



(13) 



The sum and the integral in (13) largely compensate. 
This may be seen by formally introducing an upper 
boundary fcmax = ^max, which makes the integral and 
the sum finite but drops out from their difference. The 
main contribution to this difference comes from the inte- 
gral at small fc < fcs, where ks is the radius of the Bril- 
louin zone approximated by a circle. One has kg = Ann 
for all vortex lattice symmetries. 

For A larger than the vortex spacing a, one has 
AK ^ 1 for all K. More precisely, the condition is 
= 16TT^/{3a^) > A-^ thus 27rA > a is sufficient. 
In this case the general result (13) may be evaluated ap- 
proximately by keeping just the integral over < fc < fc^ 
and noting that the remaining integral for fc > fc^ is com- 
pensated by the sum over K ^ 0. The final result for the 
curvature in this limit is then 



K','oo(o) = 



2fiQA 2fioX^ 



(14) 



IV. RESULTS FROM ELASTICITY 



The spring constant (14) may be calculated also from 
the elasticity theory of the vortex latticer^i^i^ see Ap- 
pendix A. In this case one requires nonlocal elasticity 
with a dispersive compressional modulus cii(fc). The re- 
sult coincides with Eq. (14). 

In Appendix A also is calculated the spring constant 
K;'iast (0) felt by the central vortex when the surrounding 
vortices are not pinned but are held in place only by the 
interaction with their neighbors. In this case, the elastic 
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response of the central vortex depends only on the shear 
modulus C66 of the triangular vortex latticer^i^ which is 
not dispersive. Thus, local elasticity theory is sufficient 
for this problem. However, now an inner cut-off radius 
fcmin — ^/ Rd IS required for the k integral, which is taken 
from the finite radius Rd of the thin disk containing the 
vortex lattice. It is assumed that the vortices arc pinned 
at the edge of the disk, at r = Rd- The resulting elastic 
spring constant is then 



V" 

elast 



(0) 



K%(0) 

VoAln(i?dfcs) 21n(7ri?rf/a) 



(15) 



This V^'iast('^) smaller than the spring constant 
K"oo(0) of the rigidly pinned vortex lattice by a factor 
0.5/ ln(7r i?d/a). A constant force /o on the central vor- 
tex thus displaces this vortex a distance mq that increases 
by a factor 21n(7ri?d/a) when pinning of the surrounding 
vortices is switched off. 

In all the above expressions the vortex core radius 
Tc ~ ^ and the coherence length f were assumed to 
be smaller than any other relevant length, in particular 
^ <^ h. — \^ / d. This is the London limit. Finite ^ is easily 
considered by taking the k integral not to infinite k but 
to some upper cut-off value femax ~ This cut-off re- 
moves the logarithmic divergence of V{r) (x ln(2.27A/r), 
Eq. (4), and smoothes it over the core radius ^- This 
smoothing may also be performed by repla cing in the log- 
arithm r by (r^ + r^)^/^ as shown in Ref. [slsMll. The 
exact numerical solution of the GL equations for a vortex 
and the periodic vortex lattice in films of finite thickness 
is given in Ref. d. The London solution for a vortex in 
films of finite thickness is presented e.g. in Ref. 0. 



V. VORTEX LATTICE OF FINITE SIZE 

Expressions (12) to (14) apply to an infinitely extended 
vortex lattice. The finite size of the superconducting film 
can be approximately accounted for by taking the sums 
X^R^^o ™ -^^'-l^- i^'^) over a finite area only, e.g., over a 
circular area of radius Rd — Na (disk radius), |R| < Rd, 
containing « ttN'^ vortices. The approximate potential 
that a central vortex sees in the presence of a regular 
lattice of pinned vortices in such a thin disk is thus 



T4(r,7V) =^F(r-R), 0<|R|<iVa. 



(16) 



R 



For TV — » oo, Eq. (16) coincides with Eq. (12), and the 
curvatures coincide, V;"(0,oo) = V^'{0), Eq. (13). For 
finite TV, the interaction and its curvature are reduced. 
To visualize this we define the dimensionless curvature 



f{a/A,N)^V^'iO,N)/V^'(0), 



(17) 



with /(0,oo) = 1. This function is shown in Fig. 2. 
One can see that the dense-lattice limit, Eq. (14), is a 
good approximation when a < 2A and N > 20. For 




FIG. 2: The reduced curvature fix,N) = V^'{0,N) ■ 
(2/xoA/$oB), Eq. (17), of the potential exerted on a central 
Pearl vortex by the surrounding nN^a^ pinned vortices within 
a circular area of radius Na. The function / is plotted ver- 
sus the reduced vortex spacing x — a/A for several values of 
Af = 2 to iV = oo. While f{x, N) (upper plot) for finite N has 
a maximum and vanishes at small and large x, the product 
f{x, N) ■ X for all A'' is monotonically increasing with x. 



larger a and smaller A'' the reduced curvature /(a/A, N) 
decreases. At small values of a; = a/A < 2 ior N > 2, 
f{x) has a maximum and decreases to zero when a; — + 
because then also the disk radius Rd — aN goes to zero. 



VI. VORTEX INTERACTION IN FILMS OF 
FINITE SIZE 

The interaction potential and force between two vor- 
tices in thin films of finite size and arbitrary shape de- 
pends on the size and shape of the film. Moreover, it 
depends on the positions of both vortices, V = ^(ri, r2), 
while for infinite films (Sec. I) V{r) depends only on the 
distance r = |ri — r2|. These dependences are obvious 
from the fact that the interaction potential between two 
vortices has to vanish when one vortex approaches the 
edge of the film, where its magnetic flux goes to zero since 
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FIG. 3: Current stream lines of a pair of vortices, coincid- 
ing with the contour hues of the vortex-vortex interaction. 
Shown are the contours of the logarithm of the interaction 
potential, Eq. (18), for a square thin superconducting film 
of size 2w X 2w with slit and central square hole, containing 
a pair of Pearl vortices (due to the mirror symmetry of this 
computation) at Xi = 0.54«;, yi = ±Xi; 88 x 88 grid points are 
used and 20 contours. The force exerted by this vortex pair 
on a third probing vortex at position x, y acts perpendicular 
to these contour lines. Top: 2D penetration depth A = 
(ideal screening). Bottom: A = 0.15u;. 



supercurrents cannot circulate around a vortex core po- 
sitioned on the edge. Furthermore, the interaction force 
in general is not a central force but close to the film edge 
the force on the second vortex is directed perpendicular 
to the edge no matter where the first vortex is positioned. 
This is so since the current generated by a vortex fiows 
along the film edge close to the edge, see Fig. 3. Though 
the vortex interaction energy vanishes at the edges, the 
interaction force has a peak there which increases with 
decreasing A. Only when the two vortices are far from 
all film edges and are close to each other, then the cen- 




FIG. 4: 3D plot of the interaction of a probing Pearl vortex 
positioned at {x, y) with a pair of vortices sitting at Xi = 0.4w, 
yi — ztO.Sw in a thin square film of size 2w x 2w with slit and 
hole. Plotted is the interaction potential, Eq. (18), in units 
^o/itJ'O'w), versus x,y in units of w; 46 x 46 grid points are 
used. Top: A = (ideal screening). Bottom: A — 0.15w. 



tral Pearl interaction potential V{r) of Sec. I is a good 
approximation, see Sec. VII. 

The vortex interaction in realistic films of finite ex- 
tension can be computed by the method described in 
Ref.Hilli. One has V{r„rj) = V^j = Vj, = -$o.g.(rj) 
where gi{r) is the stream function of the 2D sheet cur- 
rent density Ji(r) caused by a vortex centered at r^. 
In general one has J(r) — — z x Vg(r) = V x (zg) — 
{dg/dy, —dg/dx). The function g(r) — g{x,y) has sev- 
eral useful properties listed in Sec. 2A of Ref. [2^. In 
particular, its contour lines are the stream lines of the 
sheet current, and it may be put g = on the outer edge 
of the film, which coincides with a stream line if there is 
no current fed in by contacts. 

The function g{x,y) and its sheet current 3{x,y) can 
be caused by an applied magnetic field, or by the fiux 
trapped in a hole in the film, or by vortices, and also by 
applied currents, which we do not consider here. Within 
London theory all these contributions superimpose lin- 
early. In Ref. [2^ it is shown how these sheet currents 
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FIG. 5: Same geometry as Fig. 4, but Pearl potential, 
Eq. (19). This approximate interaction potential between 
a probing Pearl vortex at {x, y) with a pair of vortices at 
Xi = OAw, Ui = rhO.Sui does not vanish at the edges of the 
square film and of its slit and hole, while the correct numerical 
interaction (18) depicted in Fig. 4 does vanish there. 

can be computed by introducing a grid of M equidistant 
or nonequidistant points and then inverting a M x M 
matrix. The inverted matrix i^A, Eq. (11) of Ref.'l^, (di- 
mension m) is closely related to the vortex interaction. 
One has 

V{r,,r,) - ^(r„r,) = -^-^ <f^ A^^M , (18) 

where Wi are the weights of the grid points (dimension 
m^). This potential is repulsive (positive) and sharply 
peaked at = Vj, and like the function 17(2;, y) it vanishes 
at the outer edge of the film. This vanishing is linear for 
A > 0, while for A = (ideal screening) g and V go to 
zero like \/6 where S is the distance to the edge. 

If the film contains one or more holes (is multiply con- 
nected) the g{x,y) can have a non-zero value inside the 
hole, g{x, y) — I = const, where / is the total current that 
circles around the hole when the hole contains trapped 
magnetic flux. If the hole is connected with the outer 
edge by a slit, then its border is part of the outer edge 
and one has 5 = along the entire (inner and outer) edge 



of the film. If the slit is short-cut, one can trap flux in the 
slit and hole. If the slit is bridged by a superconducting 
weak link, this film can be used as a Superconducting 
Quantum Interference Device (SQUID). ^^■■^^ 

Figure 3 shows the contour lines of this vortex interac- 
tion, coinciding with the stream lines of the sheet current 
3{x,y) since V (x g. Depicted are the contours of the log- 
arithm of the interaction potential, Eq. (18), between a 
probing vortex sitting at (x, y) and a pair of vortices at 
(xijibj/i) since our computation assumes mirror symme- 
try about the x axis, \n\V(xi,yi\ x, y) -\- V(xi, —yi\ x, y)]. 
In the regions where V is very small, this logarithm shows 
more contours than a linear contour plot would show. 

Figure 4 shows the same interaction depicted as a 3D 
plot. Both Fig. 3 and Fig. 4 show the geometry of a thin 
superconducting square of size 2w x 2w with a central 
square hole and an open slit, which may be used as a 
SQUID. Depicted are the cases A = (ideal screening) 
and A/w = 0.15. One can see that for A = the potential 

V and the stream function g at the edges go to zero with 
vertical slope, V (x g cc \/S {S — distance to the edge), 
and for flnite A, V oc g oc S vanish linearly at the inner 
and outer edges of the fllm. In Fig. 4, a grid of only 46 x 46 
nonequidistant points is chosen such that no grid points 
sit directly on the fllm edges. Therefore, on the grid 
points closest to the outer edges of the film the depicted 

V is not exactly zero even when it is zero on the very 
edges. 

Figure 5 shows, for the same geometry and grid as 
in Fig. 4, the approximate vortex-vortex interaction ob- 
tained from the Pearl potential of Sec. 1 valid for an 
infinite film. To allow comparison with the correct nu- 
merical potential of Fig. 4, which has mirror symmetry, 
we symmetricize this approximate potential by plotting 

V{r,r,) = V[{{x - x,)^ + {y - y,f + e)'^^] 

+ V[{{x-x,)^ + {y + y,f + ef^] , (19) 

where V{r) is the Pearl interaction of Sec. I, e.g., Eq. (6). 
The cut-off e = 4 • 10~'^w'^ was chosen to reach the same 
peak heights as in Fig. 4. Note that the numerical in- 
teraction potential, Eq. (18) and Fig. 4, has a natural 
cut-off and finite peak heights, which are related to the 
spacing of the grid points. The peaks of the numerical V 
(Fig. 4) are well approximated by the cut-off Pearl peaks 
of Eq. (19) (Fig. 5) if they are far from all edges. 

As the peak position r.^ of V{r, r^) approaches the film 
edge, the amplitude of the correct potential and its peak 
height decrease and finally vanish as the edge is reached. 
This means that the interaction energy of two vortices is 
strongly reduced when both vortices are close to the film 
edge. The interaction force, however, i.e., the slope of V, 
may still be large, especially when A is small and r or 
or both are close to the edge. 

The VS behavior of the vortex interaction near the film 
edge can be seen in Fig. 5 of Ref. H^, which shows this 
potential (there called integral kernel K) for an infinite 
(along y) thin strip with A = 0, containing a dense row 
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of vortices (same Xi, many yi); also shown there is the 
product K{x, Xi) ■ \x — Xi\^ that enlarges the y/S shape at 
the edges. This idcal-scrccning potential has a In |a; — x^j 
singularity obtained by integrating the l/|r — r^l singu- 
larities of the Pearl potential (for A = 0) along the posi- 
tions Ui- 



A = 




FIG. 6: The ratio of the exact numerical interaction po- 
tential between two vortices and the Pearl interaction, 
Vnuni{x,y)/Vpciiri{^,y), for a thin square of size 2™ x 2w with 
one vortex positioned at the square center (0, 0) and the other 
vortex at (x,y). Top: A = 0. Bottom: A = 0.3w. 



VII. RATIO OF NUMERICAL AND PEARL 
POTENTIALS 

The analytical Pearl potential (19) (Fig. 5) does not 
depend on the shape of the film, in particular, it does 
not vanish at the outer and inner edges as it should. If 
an approximate analytic interaction potential is needed, 
e.g., for computer simulation of vortex motion and vor- 
tex pinning in films, one may construct this for simple 
film shapes like circular disks or rectangles or long strips, 
by multiplying the Pearl potential by a factor f{x,y) 
that vanishes at the film edges with a A dependent slope 
and that looks like a conical mountain whose maximum 



height for A = reaches unity and for A > is somewhat 
lower than unity. This numerical finding means that for 
A > the correct, shape-dependent interaction is smaller 
than the ideal Pearl potential. The physical reason for 
this is that for finite films the magnetic field lines of a vor- 
tex do not extend to infinity but return around the film 
edges, thereby reducing the sheet current that causes the 
vortex vortex interaction, sec discussion before Eq. (1). 

Two examples for this factor f{x,y) = 
^um(a;, y)/Vpeari(a;, y) , i.e., the ratio of the exact 
numerical potential and the ideal Pearl potential, arc 
depicted in Fig. 6 for a superconducting thin film square 
with the first vortex positioned at the center (0, 0) and 
the second vortex at (x, y). This ratio has approximately 
conical shape, i.e., the correction factor Kium/l^eari 
nearly linearly goes from unity to zero as the second 
vortex moves from the center to the edge of the square. 
It can be nicely seen in Fig. 6 that near the edges 
Kium(a;,y) vanishes cx 5^^'^ when A = and oc 5 when 
A > 0. These plots also show the somewhat unexpected 
result that the deviation of the exact interaction poten- 
tial from the Pearl potential can be considerable even 
when both vortices are far from the film edges. 



VIII. SUMMARY 

It is shown that the interaction potential between two 
Pearl vortices in a superconducting thin infinite film can 
be well approximated by a simple logarithm at all dis- 
tances, Eq. (6) and Fig. 1. A general argument is given 
why this interaction at large distances r » A has the uni- 
versal behavior V{r) = $o/(Mot), independent of the 
2D magnetic penetration depth A = X^/d. Explicit ex- 
pressions are given for the interaction of one central Pearl 
vortex with an infinite regular lattice of pinned Pearl vor- 
tices or with the image vortices in numerical simulations 
that use periodic boundary conditions. The curvature 
of this interaction (spring constant of the central vortex) 
is presented both for this regular rigid lattice, and for 
the elastically deformable lattice of Pearl vortices. The 
modification of these results by the finite extension of 
the vortex lattice in finite films (e.g., disks) is given, see 
Fig. 2. Finally, it is shown how the interaction of the vor- 
tices in thin films of finite size and arbitrary shape can 
be computed. This general vortex-vortex interaction po- 
tential depends not only on the distance r = |ri — r2| 
of the two Pearl-like vortices but on both positions ri 
and r2. The interaction in general is not central; e.g., 
the interaction force on a vortex near the film edge acts 
perpendicular to this edge no matter where the other 
vortex is positioned. Anyway, this non-central interac- 
tion is caused by the other vortex and it vanishes when 
one of the two vortices approaches the edge. Examples 
for this correct interaction are shown in Figs. 3 and 4. 
Figure 5 shows the corresponding radial-symmetric and 
film-shape independent Pearl potential that, away from 
the film edges, exhibits the approximately correct peaks 
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but does not vanish at the (inner and outer) film edges 
as the correct potential does. The correct numerical vor- 
tex interaction in finite-size films everywhere in the film 
is weaker than the ideal, infinite-film vortex interaction 
potential of Pearl, except when both vortices are close to 
each other near the middle of the film. This is shown in 
Fig. 6 for a square film with one vortex at the center. 
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APPENDIX A: ELASTICITY THEORY 

The linear elastic restoring force experienced by a dis- 
placed vortex in a 2D lattice of pinned or pin-free vortices 
may also be calculated from the theory of elasticity of the 
vortex lattice. Within this 2D problem the vortex lattice 
is described by its uniaxial compressional modulus cn 
and by its shear modulus cgg, which is typically much 
smaller, cqq <^ cn. Within London theory one has for 
the triangular lattice of parallel Abrikosov vortices22i^i^ 



Cii(fc) = 



1 + /c2A2 ' 



C66 



167r/ioA2 



(Al) 



where k = [k^^ky) is the wave vector of the periodic 
displacement field. The k dependence (dispersion) of 
Ciiik) means that the elasticity is nonlocal. One has 
the ratio cii(0)/c66 = SfeK^ ^ (327r/\/3)(A/a)2 > 1, with 
h = B/Bc2 and Bc2 — *i'o/(27r^^) the upper critical field 
of the superconductor. Expressions (Al) are valid when 
the magnetic fields of the vortices overlap, i.e., for vor- 
tex spacings a <^ 2'kX equivalent to 1/{2k'^) <C 6 ^ 1. 
At very small inductions, approximate expressions may 
be obtained by considering only nearest neighbor inter- 
action, yielding cn (x cee oc exp(— a/A). 

Interestingly, for the 2D lattice of Pearl vortices in thin 
films the same moduli cn, cge, Eq. (Al), apply (referred 
to unit volume). Namely, in the limit of thin films with 
thickness d <C A in the London limit k ^ 1 one has^ 
for the moduli per unit area dcii{k) = dB^ / {^ok^ X^) = 
SV(Mofc^A) (for k^X^ » 1) and dcee = B$o/(167r^oA). 
Thus, the below results apply both to the lattice of par- 
allel Abrikosov vortices with length d ^ a and to thin 
films with thickness d ^ A. 

We introduce the vortex displacements = u(Rj^) 
and their Fourier transforms u(k) — [uj;(k), Uy(k)], 



u.= / ^u(k)e»''^ u(k)=5:^, 
Jbz 47r2 ^ n 



ikR, 



(A2) 



Here n = B/^q is the vortex density and the k integral is 
over the first Brillouin zone (BZ) with area An'^n — irk^, 
where fcs = (ATm)^^^ w 7r/a is the radius of the circle 



that may be used to approximate the hexagonal BZ. At 
the maximum k — ks one has cii(fcB)/c66 = 4. The 
elastic energy per unit area (hence the factor d) is 



d f d^k 



(A3) 



In it ^^^(k) with a, (3 — x,y is the elastic matrix, 
$a/3(k) = (cu - ce6)kak/3 + Sapcuak"^ with k'^ = kl + k^. 
Explicitly one has 



^xx = ciikl + caeky, ^yy = caekl + cnky 

^xy = ^yx = (Cll - Cf,Q)kj:ky , 



(A4) 



the determinant of <I>Q^(k) is -D(k) — k'^cuCee, and 



the inverse matrix has the elements ($ 



($- 



/yy 



<P^ JD, and ($- 



Jxy 



($- 



/yx 



^xylD. 



We shall now apply this formalism to two problems. 

First, we consider the problem of Eq. (12): all vortices 
are rigidly pinned but the free vortex at the origin is 
shifted by a distance uq along x, thus — uoS,yfi'k. 
Inserting this into (A2) yields u(k) — (^io/?^)x and the 
elastic energy (A3) becomes 

Since cn ^ cqq one has from Eqs. (Al) and (A4) ^xx ~ 
cii(fc)fc2 « (BVAio)fc^/(l + fc'A2) « (SVMoA2)fc2/A:2. 
Averaging k'^/k'^ — cos'^ ip over the angle ip yields i. 
With the BZ area Air'^n we thus obtain the potential cur- 
vature 



K'.'oo(0) 



B^d 



BZ 



d2fc fc| 
4^ fc2 



B<^od 
2moA2 



(A6) 



This coincides with Eq. (14). 

Second, we consider the case when all vortices in a 
cylinder or disk with radius Rd can move freely, only the 
vortices at the surface or edge r = Rd are pinned to en- 
sure zero displacement there. We calculate the elastic 
force /o that is required to move the central vortex a 
distance uq, thereby deforming the vortex lattice. Thus, 
the applied forces on the vortices are all zero except for 
the force acting on the central vortex, fj/ = foS^fi'X- In 
general, we define the forces per unit length (the to- 
tal force on the i/th vortex of length d is di,y) and their 
Fourier transform f(k) = [/a;(k), (k)] by 



BZ 



fk 
4^ 



Ell/ 



(A7) 



The elastic energy (A3) may also be written as 
d X - 

-felast = 2 / y 

V 



d f d'k 



'BZ 



fk 

4^ 



/„(k)$;i(k)/^(k). (A8) 
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Displacements and forces arc related by 

n/„(k) = $„/3(k)u;(k), 
t.„(k)=n$;^(k)/^(k). 



(A9) 



For our example = /o(5i/.ox one has f(k) = (/o/n)x 
and the elastic energy becomes, with w ^^/(fc^cee), 



. (AlO) 



For the k integral we have used a lower cut-off R^^ and 
the upper limit fcs « w/a. The curvature of the elastic 
potential, ^4'iast(0)' vortex at the origin feels is 

defined by 



1 d^fS 



2 ^JLt(O) ' 



(All) 



since the central force is dfo — ^^,"3^3^(0)^0, i.e., K'iast(O) 
is a spring constant. Comparing Eqs. (AlO) and (All) 
and using the shear modulus (Al) we obtain 



K.ust(0) = 



B^od 



vX.o(o) 



(A12) 



Thus, when the surrounding vortices are allowed to re- 
lax elastically, the displacement of the central vortex, on 
which a force acts, is increased by a factor 2ln{nRd/a). 



^ J. Pearl, Appl. Phys. Lett. 5, 65 (1964). 

2 J.-C. Wei and T.-J. Yang, Jpn. J. Appl. Phys. 35, 5696 

(1996). 

^ D. Yu. Irz, V. N. Ryzhov, and E. E. Tareyeva, 
Phys. Lett. A 207, 374 (1995). 

G. Carneiro and E. H. Brandt, Phys. Rev. B 61, 6370 
(2000). 

^ E. H. Brandt, Phys. Rev. B 71, 014521 (2005). 

® P. S. Deo, V. A. Schwcigcrt, F. M. Pccters, and A. K. 

Gcim, Phys. Rev. Lett. 79, 4653 (1997). 

V. A. Schwcigert and F. M. Peeters, Phys. Rev. Lett. 83, 

2409 (1999). 

* V. A. Schweigert and F. M. Peeters, Physica C 332, 266, 
426 (2000). 

^ B. J. Baelus, L. R. E. Cabral, and F. M. Peeters, Phys. 

Rev. B 69, 064506 (2004). 
^° D. S. Golubovic, M. V. Milosevic, F. M. Peeters, and V. 

V. Moshchalkov, Phys. Rev. B 71, 180502(R) (2005). 
" V. R. Misko, B. Xu, and F. M. Peeters, Phys. Rev. B 76, 

024516 (2007). 

Guo-Qiao Zha, Shi-Ping Zhou, and Bao-He Zhu, Phys. 
Rev. B 73. 092512 (2006). 

B. J. Baelus and F. M. Peeters, Phys. Rev. B 65, 104515 

(2002). 

G. R. Bcrdiyorov, B. J. Baelus, M. V. Milosevic, and 

F. M. Peeters, Phys. Rev. B 68, 174521 (2003). 

G. R. Berdiyorov, M. V. Milosevic, and F. M. Peeters, J. 
Low Temp. Phys. 139, 229 (2005). 

B. J. Baelus, A. Kanda, N. Shimizu, K. Tadano, Y. 
Ootuka, K. Kadowaki, and F. M. Peeters, Phys. Rev. B 
73, 024514 (2006). 
" A. V. Silhanek, W. Gillijns, A. Volodin, V. V. 
Moshchalkov, M. V. Milosevic, and F. M. Peeters, Phys. 
Rev. B 76, 100502 (R) (2007). 



M. V. Milosevic and F. M. Peeters, Phys. Rev. Lett. 93, 
267006 (2004); Phys. Rev. B 68, 024509 (2003). 
G. R. Berdiyorov, B. J. Baelus, M. V. Milosevic, and F. 
M. Peeters, Phys. Rev. B 68, 174521 (2003). 
^° G. R. Berdiyorov, M. V. Milosevic, and F. M. Peeters, J. 
Low Temp. Phys. 139, 229 (2005). 

M. V. Milosevic, S. V. Yampolskii, and F. M. Peeters, 
Phys. Rev. B 66, 024515 (2002). 
22 L. R. E. Cabral and F. M. Peeters, Phys. Rev. B 70, 214522 

(2004) . 

A. D. Hernandez and A. Lopez, Phys. Rev. B 77, 144506 
(2008). 

E. SardcUa, P. Noronha Lisboa Filho, and A. L. Malvezzi, 
Phys. Rev. B 77, 104508 (2008). 
E. H. Brandt, Phys. Rev. B 64, 024505 (2001). 
E. H. Brandt, Phys. Rev. B 72, 024529 (2005). 
M. Abramowitz and I. A. Stegun, eds.. Handbook of Math- 
ematical Functions (National Bureau of Standards, Wash- 
ington, 1967). 

^® J. R. Clem, J. Superconductivity 17, 613 (2004). Note that 

Clem uses a different definition of A = 2A^/d. 

E. H. Brandt, phys. stat.sol.(b) 77, 551 (1976). 
^° E. H. Brandt, J. Low Temp. Phys. 26, 735, 709 (1977). 

R. Labusch, phys. stat.sol. 32, 436 (1969). 
^2 J. R. Clem, J. Low Temp. Phys. 18, 427 (1975). 
3^ E. H. Brandt, Phys. Rev. Lett. 78, 2208 (1997). 

A. Yaouanc, P. Dalmas de Reotier, and E. H. Brandt, Phys. 

Rev. B 55, 11107 (1997). 
35 J. R. Clem and E. H. Brandt, Phys. Rev. B 72, 174511 

(2005) . 

3® E. H. Brandt, J. de Physique 48, Colloque C8, 31-43 
(1987). 



